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Through-plane  liquid  water  distributions  recently  visualized  by  Manahan  et  al.  [1]  and  Turhan  et  al.  [2] 
using  the  neutron  radiography  (NR)  technique  show  that  the  peaks  of  the  water  distributions  are  located 
near  the  center  of  a  gas  diffusion  layer  (GDL).  We  suggest  that  the  distinctive  water  profiles  are  caused  by 
incomplete  polytetrafluoroethylene  (PTFE)  treatment  of  the  GDL  and  the  resultant  spatial  variation  of 
GDL  wettability  in  the  through-plane  direction.  Based  on  this  hypothesis,  we  improve  the  macroscopic 
two-phase  fuel  cell  model  to  describe  two-phase  transport  through  GDLs  with  variation  of  spatial 
wettability  [3].  The  proposed  model  successfully  reproduces  the  shape  of  through-plane  water  profiles 
obtained  from  the  NR  experiments  [1,2].  Therefore,  the  centrally  located  liquid  saturation  peak  in  the  GDL 
can  be  attributed  to  incomplete  PTFE  treatment  of  the  GDL.  This  occurs  because  liquid  water  is  more 
easily  accumulated  in  the  relatively  hydrophilic  GDL  pores  encountered  in  the  inner  GDL  region  (rather 
than  the  outer  GDL  region)  due  to  its  incomplete  PTFE  treatment.  Our  results  indicate  that  the  overall 
characteristics  of  liquid  water  distribution  in  a  GDL  under  an  inhomogeneous  wetting  condition  can  be 
macroscopically  predicted  using  the  two-phase  model  presented  here. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  gas  diffusion  layer  (GDL),  usually  made  of  carbon  paper  or 
carbon  cloth,  plays  a  crucial  role  in  the  water  management  of 
polymer  electrolyte  fuel  cells  (PEFCs).  Because  the  carbon  paper 
and  carbon  cloth  GDL  materials  inherently  exhibit  hydrophilic 
behavior,  they  are  treated  with  wet-proofing  agents  such  as  poly¬ 
tetrafluoroethylene  (PTFE)  to  more  efficiently  remove  product 
water  during  PEFC  operation.  However,  after  the  wet-proofing 
treatment,  the  PTFE  loading  may  not  be  uniform  throughout  the 
GDL  due  to  the  complex  pore  structure  of  the  GDL  [4].  Furthermore, 
a  gradual  loss  of  GDL  hydrophobicity  due  to  GDL  aging  also  exists  as 
reported  in  the  literature  [5,6]. 

Recently,  Manahan  et  al.  [1]  and  Turhan  et  al.  [2]  used  the 
neutron  radiography  (NR)  technique  to  visualize  through-plane 
water  profiles  along  the  thickness  of  a  GDL.  As  shown  in  Fig.  1, 
their  NR  data  show  that  the  peaks  of  the  water  distributions  are  all 
located  near  the  center  of  the  GDL  thickness.  Assuming  uniform 
GDL  properties,  they  concluded  that  this  trend  resulted  from 
a  combined  effect  of  capillary  flow  and  phase  change-induced  (PCI) 
flow.  However,  we  emphasize  that  this  type  of  water  profile  can  be 
obtainable  solely  by  capillary  action,  as  long  as  porous  properties 
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such  as  porosity  and  contact  angle  are  spatially  non-uniform  inside 
the  GDL.  Pore-scale  simulation  results  described  in  the  literature 
demonstrate  that  the  spatial  wettability  variation  throughout 
a  GDL  significantly  alters  capillary  transport  and  the  resultant  water 
profile  inside  the  GDL  [7]. 

Previously,  we  developed  a  two-phase  fuel  cell  model  in  which 
a  new  feature  regarding  the  effect  of  the  spatial  variation  of  the  GDL 
porous  properties  was  modeled  and  implemented  [3].  From  our 
simulation  results,  we  found  that  the  wettability  variation  through 
the  GDL  thickness  considerably  altered  the  water  transport  char¬ 
acteristics  and  resultant  water  distribution  in  the  GDL.  The  present 
study  is  an  extension  of  our  previous  work  [3],  with  the  goal  of 
reproducing  the  shape  of  through-plane  water  profiles  obtained  by 
Manahan  et  al.  [  1  ]  and  Turhan  et  al.  [2  ].  By  reasonably  assuming  the 
wetting  property  (contact  angle)  variation  along  the  GDL  thickness, 
our  simulation  results  clearly  demonstrate  that  a  water  peak  can 
exist  near  the  center  of  the  GDL  solely  driven  by  capillary  transport. 


2.  Numerical  model 

2.1.  Two-phase  fuel  cell  model 

The  one-dimensional  (1-D)  two-phase  model  presented  in  our 
previous  study  [3]  was  used  to  predict  the  experimental  liquid 
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Fig.  1.  NR  data  for  through-plane  liquid  water  distributions  in  the  anode  and  cathode  GDLs.  Data  were  replotted  from  (a)  Manahan  et  al.  [1]  (SGL  10BB  type  GDL  with  250  pm 
thickness;  active  area,  Amem  =  5  cm2;  cell  temperature,  Tceii  =  65  °C;  anode  and  cathode  inlet  humidification,  RHa/RHc  =  100%;  anode  and  cathode  stoichiometric  ratios,  £a/£c  =  2.0/ 
2.0)  and  (b)  Turhan  et  al.  [2]  (SGL  10BB  type  GDL  with  300  pm  thickness;  active  area,  Amem  =  17.2  cm2;  anode  and  cathode  inlet  pressures,  Pa/Pc  =  1  atm;  anode  and  cathode  inlet 
humidification,  RHa/RHc  =  100%;  anode  and  cathode  stoichiometric  ratios,  fa/£c  =  1.5/2.0). 


water  profile  through  the  GDL  thickness.  The  model  is  based  on 
a  multiphase  mixture  (M2)  approach,  and  considers  the  effects  of 
immobile  liquid  saturation  and  the  spatial  variation  of  GDL 
wettability  on  the  liquid  water  distribution  inside  a  GDL.  The 
computational  domain  for  the  model  includes  the  anode  GDL,  the 
catalyst-coated  membrane  (CCM),  and  the  cathode  GDL.  Therefore, 
gas  channels  (GCs)  and  bipolar  plates  (BPs)  for  both  the  anode  and 
cathode  sides  are  excluded  to  simplify  numerical  analysis. 
Regarding  the  water  transport  across  the  membrane,  all  three 
modes  comprising  the  hydraulic  pressure  gradient,  the  membrane 


water  content  gradient,  and  the  electro-osmotic  drag  due  to  proton 
flux  were  rigorously  taken  into  account  for  a  more  accurate 
prediction  of  water  transport  and  accumulation  inside  a  cell.  The 
computational  domain  and  the  transport  processes  for  the  present 
model  are  schematically  shown  in  Fig.  2. 

2.2.  Governing  and  constitutive  equations 

The  present  two-phase  model  is  governed  by  conservation  of 
mass,  momentum,  and  species.  All  the  governing  equations  of  this 
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Fig.  2.  Schematic  of  computational  domain  (red  dashed  line)  and  electrochemical  reactions/transport  phenomena  considered  in  the  model.  HOR  and  ORR  denote  hydrogen 
oxidation  reaction  in  anode  side  and  oxygen  reduction  reaction  in  the  cathode  side,  respectively  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.). 


model  are  listed  in  Table  1.  Additionally,  the  transport  properties  in 
the  membrane  and  the  cell  design/material/physical  properties  are 
summarized  in  Tables  2  and  3,  respectively.  Readers  are  referred  to 
Ju  [3]  for  more  detailed  elaboration  of  the  model.  Since  the  focus  of 
the  present  study  is  on  water  transport  through  GDLs  with  varia¬ 
tion  of  spatial  wettability,  the  water  transport  equation  for  this  type 
of  GDL  is  described  here  in  more  detail. 

The  steady  state  1-D  water  conservation  equation  in  the 
through-plane  direction  (x)  can  be  written  as  [3] 


d_ 

dx 


(YwPmwu) 


d_ 

dx 


[(m^-m^)/] 


(12) 


It  should  be  noted  that  in  the  M2  model,  a  diffusive  mass  flux  of 
each  phase  relative  to  the  whole  multiphase  mixture  is  defined,  and 
j1  in  the  second  term  on  the  right-hand  side  of  Eq.  (12)  denotes  the 
diffusive  mass  flux  of  the  liquid  phase  relative  to  the  whole 
multiphase  mixture  as  follows  [11]: 


jl  =  plul  -  Xlpu 


(13) 


In  Eq.  (14),  the  capillary  pressure  Pc,  which  is  defined  as  the 
difference  between  the  gas  and  liquid  pressures,  is  expressed  as  [3] 


Pc  =  ps  _  pi  =  (TCOS  e(£)  1/2J(sr)  (15) 

where  e}  I<,  and  0  are  the  porosity,  permeability,  and  contact  angle  of 
the  given  porous  GDL,  respectively.  The  Leverett  function,  J(sr), 
which  is  a  dimensionless  capillary  pressure,  can  be  expressed  as 
a  function  of  the  liquid  saturation  sr  for  the  given  hydrophobic  GDL 
as  follows: 

fo  ifSr<Sim  nfis 

\  1.417sr -2.120s2  +  1.263s?  ifsr>sim  (  > 

In  Eq.  (16),  S[m  is  the  immobile  liquid  saturation  and  sr  is  the 
liquid  saturation  representing  the  volume  fraction  of  GDL  pores 
occupied  by  transportable  liquid  as  follows: 


Sr  = 


V'/Vpore-Si 
1  Sim 


(17) 


Therefore,  the  second  term  on  the  right-hand  side  of  Eq.  (12)  (i.e. 
the  capillary  transport  term)  is  due  to  the  relative  motion  between 
the  gas  and  liquid  phases  and  is  proportional  to  the  capillary 
pressure  gradient  as  follows: 

(mw  -  =  (1  -  mh)  (14) 


It  should  be  emphasized  here  that  the  standard  Leverett 
J-function  form  in  Eq.  (16)  was  formulated  to  characterize  the  liquid 
transport  in  more  uniform  porous  media  such  as  soils  and  rocks 
[12,13].  A  more  realistic  capillary  pressure-saturation  study 
focusing  on  real  fuel  cell  GDL  was  conducted  by  Gostick  et  al.  [14] 
and  Kumbur  et  al.  [15,16].  In  particular,  Kumbur  et  al.  [15,16]  con¬ 
ducted  capillary  pressure-saturation  measurements  under  typical 
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Table  1 

Governing  equations  of  the  1-D  two-phase  fuel  cell  model. 


Mass 


Momentum 

Species 


Governing  equations 


For  the  anode  GDL: 


/  , 


I 


pu  =  2pMH2  +  |ndp  -  Dj 
For  the  cathode  GDL: 

/  „  I  „ 


SI  1X161X1  / 
_ ( 

dA\ 

EW  \ 

to) 

'  1 

\  I 

nd  +  2 

/Cmem/dpl 

Mw - j —  “jtt" 


dx 


-  D17 


EW 


Mw 


For  the  anode  and  cathode  GDLs: 
KdP 
v  dx 


Kmem  / dpl 

jA  Id* 


pu 


Water  transport  in  the  anode  GDL: 

7wpm„u  =  PsD%fs^  (mi)  +  (mi  -  m*,)  /  +  ["dp  -  Dwem^- Mw  - 

Water  transport  in  the  cathode  GDL: 

jwpmwu  =  Pg£>weff^  (ml)  +  (ml  -  m\ v)  j1  +  [(nd  +  ^  -  Dwem^w“(^) 


Mw  - 


JCmem  / dP1 
vl  l  dx 


Water  transport  across  the  membrane: 


„dA 


Mw-- 


ndp  I  Ml 


d  /omei 

~  (  r _ nmem _  _ 

dx  V  EW  w  dx/  iVlw  dx 
Oxygen  transport  in  cathode  GDL: 

10lpm0lu  =  ^ofg^Hh)  +m§/ 


d  /KmemdP>\ 

'+dx^  j/>  dx  J  ~ 

->o2 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
(7) 


fuel  cell  compressed  GDL  and  operating  temperatures  (20,  50,  and 
80  °C).  They  compared  their  experimental  measurements  and 
capillary  pressure-saturation  correlation  test  data  obtained  with 
the  traditional  Leverett  function.  According  to  their  comparison, 
the  deviation  between  the  Leverett  function  and  new  correlations 
is  not  significant,  especially  for  the  practical  range  of  liquid  satu¬ 
ration  (0  <  s  <  0.5)  in  GDL  during  fuel  cell  operation.  Therefore,  we 
employ  the  standard  Leverett  J-function,  i.e.  Eq.  (16),  for  this  study. 

Eq.  (14)  shows  that  capillary  pressure  Pc  is  a  basic  parameter 
characterizing  liquid  water  transport  in  a  porous  GDL.  By 
substituting  the  capillary  pressure-related  terms  in  Eqs.  (15)— (17) 
into  Eq.  (14),  the  final  form  of  the  capillary  transport  term 
described  by  Eq.  (14)  can  be  written  as  follows: 


(1  -  ”4) 


X'Xs<r(eK)1/2 

V 


cos  0 


d J  d sr 
dsr  dx 


second  term  on  the  right-hand  side  of  Eq.  (18)  is  the  distinctive 
feature  of  the  proposed  model  and  is  necessary  to  describe 
two-phase  transport  through  a  porous  medium  with  inhomoge¬ 
neous  surface  wettability.  However,  the  second  term  is  not  included 
in  most  macroscopic  two-phase  fuel  cell  models  introduced  in  the 
literature. 

2.3.  Boundary  conditions 

As  indicated  in  Fig.  2,  the  aforementioned  one-dimensional 
model  requires  boundary  conditions  for  water  and  oxygen  species 
at  the  GC/GDL  interfaces.  The  interfacial  mass  fractions  for  oxygen 
and  water  species,  iti^gdl/gg  are  determined  by  the  cell  operating 
pressure  (P),  temperature  (T),  GC  relative  humidity  (RHgc).  and 
interfacial  liquid  droplet  coverage  S[nt  due  to  channel  flooding 
conditions. 


+J(Sr) 


d(cos  0)1/2 
dx 


(18) 


where  the  first  and  second  terms  on  the  right-hand  side  of  Eq.  (18) 
represent  capillary  transport  driven  by  the  liquid  saturation 
gradient  and  the  spatial  variation  of  the  GDL  wetting  characteristic 
(contact  angle  0),  respectively.  It  should  be  emphasized  that  the 


Table  2 

Transport  properties  of  electrolyte  phase. 


Expression 
Water  activity  [8]: 

CjRuT 

psat 

Water  content: 

3  f  Ag  =  0.043  +  17.81a  -  39.85a2  +  36.0a3  for  0  <  a  <  1 
A=  1^  =  22 

Electro-osmotic  drag  coefficient  [9]: 

f  1  A  <  Ag  (a  =  1) 

d  \  2.5  A  =  A1 

Water  diffusion  coefficient  in  the  membrane  (m2  s_1)  [10]: 

r  3.1  x  10-7A(ea28;i  -  l)e(-2346/T)  for  0  <  A  <  3 
£)mem  _  I  V  J 

A  l  4.17  x  10-8A^1  +  161e-;i)e(-2346/T)  otherwise 


(B) 

(9) 

(10) 


(11) 


mw,GDL/GC  (5int  j  RHgc  ,T,P ) 


1  /  -I  \  PwMv 

P  5int  +  (1  —  Sint)  g  j 

/^5int  +  P§(1  —  Sjnt) 


P^int  +  (1 


5int) ' 


RHGC  •  Psat  •  Mw 
RuT 


Table  3 

Cell  design,  material  and  physical  parameters. 

P'Sjnt  +  P®(1  ^int) 

(19) 

Description 

Value 

Porosity  of  the  GDL,  £Gdl 

0.6 

Permeability  of  the  GDL,  XGdl 

1.0  x  10-12  m2 

Thickness  of  GDL,  5GDl 

0.3  mm 

Hydraulic  permeability  of  the  membrane,  I<mem  5.0  x  10  20  m2 

Thickness  of  membrane,  5mem 

0.018  mm 

Dry  membrane  density,  pmem 

2000  kg  rrr3 

Equivalent  weight  of  electrolyte  in  the  membrane,  EW  1.1  kg  mol  1 

Faraday  constant,  F 

96,487  C  mol-1 

Universal  gas  constant,  Ru 

8.314  J  mol  1  K  1 

Surface  tension,  a 

0.0625  N  nr1 

Liquid  water  density,  p1  (80  °C) 

972  kg  m-3 

Liquid  water  viscosity,  p1 

3.5  x  10-4  N  s  m"2 

K.  Kang  et  al.  /  Journal  of  Power  Sources  212  (2012)  93-99 


97 


m02,GDL/Gc(5inL  RH,  T,P) 


(1  Sint)Mc>2^^ 

P^int  +  PS0  —  5int) 

n  -  )M  P“Pw-PN2 
(1  Sint)Mo2 

int  +  PS0  —  5int) 


(20) 


2.4.  Numerical  procedures 

The  first  order  ordinary  differential  equations  shown  in  Table  1 
are  solved  separately  in  three  different  regions,  i.e.  anode  GDL, 
CCM,  and  cathode  GDL  20  grid  points  for  each  of  three  computa¬ 
tional  regions  was  found  to  provide  sufficient  spatial  resolution.  In 
order  to  connect  these  differential  equations  for  the  three  domains, 
an  initial  estimated  value  is  provided  for  the  water  flux  across  the 
membrane  to  set  up  the  interfacial  boundary  conditions.  An  itera¬ 
tive  procedure  is  used  to  improve  the  initial  estimate  where  the 
individually  calculated  interfacial  fluxes  in  adjoining  domains 
should  be  matched  at  the  interface.  The  iterations  proceed  until  the 
relative  error  falls  below  the  convergence  criterion  (10  7). 

3.  Results  and  discussion 

We  expected  that  the  wetting  characteristic  of  a  GDL  would  vary 
spatially  due  to  several  possible  reasons  such  as  heterogeneous  GDL 
pore  structures,  anomalies  in  the  PTFE  treatment,  surface  defects, 
impurities,  and  aging  of  the  GDL.  We  focus  on  the  effect  of  PTFE 
treatment  on  the  liquid  saturation  profile.  During  PTFE  treatment, 
a  GDL  is  soaked  in  PTFE  solution  to  enhance  its  water  removal 
capability.  As  shown  schematically  in  Fig.  3a,  it  is  possible  that  the 
percolation  of  PTFE  solution  into  a  GDL  (i.e.  mainly  driven  by 
capillary  effects)  is  not  sufficient,  causing  fewer  PTFE-coated  pores 
near  the  center  of  a  GDL  than  near  the  outer  GDL  region.  Conse¬ 
quently,  the  center  of  a  GDL  can  exhibit  less  hydrophobic  behavior 
than  the  outer  GDL  surface,  as  shown  in  Fig.  3b.  Recently,  Rofaiel 
et  al.  [4]  measured  PTFE  distributions  in  the  through-plane  direc¬ 
tion  for  paper,  felt,  and  cloth  GDLs  and  their  experimental  data 


Fig.  4.  (a)  Schematic  representation  of  spatial  variation  of  the  cathode  GDL  contact 
angle  along  the  through-plane  direction  for  Cases  1,  2,  3,  and  4.  (b)  Comparison  of 
calculated  through-plane  liquid  saturation  profiles  in  the  cathode  GDL  The  simulation 
conditions  were  at  80  °C,  1.5  A  cm-2,  and  1.5  atm  for  both  anode  and  cathode  GC  inlet 
pressures. 


exhibit  lower  PTFE  concentrations  near  the  GDL  core  region  as 
compared  to  the  outer  GDL  surface  for  the  paper  and  felt  GDLs. 

In  our  previous  study  [3],  the  effect  of  variation  in  GDL  wetta¬ 
bility  along  the  GDL  thickness  was  numerically  investigated.  Three 


Fig.  3.  (a)  Schematic  representation  describing  percolation  of  PTFE  solution  into  a  GDL  during  PTFE  treatment  process,  (b)  Possible  variation  of  GDL  wettability  along  the  through- 
plane  direction  after  PTFE  treatment. 
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cases  were  defined  by  varying  the  cosine  of  the  GDL  contact  angle 
(cos  8)  along  the  GDL  thickness.  Case  1  is  a  reference  case  repre¬ 
senting  no  variation  of  cos  6  in  the  GDL  thickness  direction  (i.e.  a 
uniform  contact  angle)  while  the  value  of  cos  8  linearly  increases 
(Case  2)  or  decreases  (Case  3)  from  the  CCM/GDL  interface  toward 
the  GC.  To  analyze  the  effect  of  incomplete  PTFE  treatment,  an 
additional  case  (Case  4)  is  newly  defined  and  simulated  in  this 
study.  As  schematically  described  in  Fig.  4a,  the  absolute  value  of 
cos  8  for  Case  4  decreases  in  the  first  half  of  the  GDL  and  then 
increases  again  in  the  second  half  of  the  GDL. 

Fig.  4b  compares  the  calculated  liquid  saturation  profile  for  Case 
4  to  those  of  Cases  1, 2,  and  3.  Since  all  four  cases  in  Fig.  4b  are  based 
on  the  high  current  density  fuel  cell  operation  (1.5  A  cm-2),  no 
liquid  droplet  coverage  at  the  cathode  GC/GDL  interface  is  assumed 
due  to  the  high  gas  flow  rate  in  GCs  and  thus  efficient  droplet 
removal  on  the  GDL  surface  [17,18].  The  liquid  saturation  peak  for 
Case  4  forms  near  the  center  of  the  GDL  because  liquid  water 
accumulation  is  more  strongly  promoted  by  the  relatively  hydro¬ 
philic  GDL  inner  region.  The  shape  of  the  through-plane  water 
distribution  predicted  by  the  proposed  model  qualitatively 
resembles  the  water  profiles  experimentally  captured  by  Manahana 
et  al.  [1  ]  and  Turhan  et  al.  [2]  using  the  neutron  imaging  technique. 

Fig.  5  shows  the  liquid  saturation  distributions  of  Case  4  at  two 
different  current  densities  of  0.2  A  cm-2  and  1.2  A  cm-2.  According 
to  the  experimental  observations  from  Yang  et  al.  [17]  and  Zhang 
et  al.  [18],  the  number  of  liquid  droplets  on  the  cathode  GDL/GC 
interfacial  surface  is  inversely  proportional  to  the  gas  flow  rate  in 
GCs  and  thus  operating  current  density.  Therefore,  we  assume  the 
values  of  interfacial  droplet  coverage  at  the  cathode  GDL/GC 
interface  for  the  cases  of  0.2  A  cm-2  and  1.2  A  cm-2  to  be  15%  and 
7.5%,  respectively.  The  liquid  saturation  peaks  observed  in  these 
cases  exists  almost  near  the  center  of  the  GDL,  which  is  driven  by 
the  assumed  contact  angle  variation  to  approximate  incomplete 
PTFE  treatment  (see  Case  4  in  Fig.  4a).  More  importantly,  these  two 
cases  exhibit  almost  the  same  level  of  liquid  saturation  peak  in  the 
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Fig.  5.  Comparison  of  calculated  through-plane  liquid  saturation  profiles  in  the 
cathode  GDL  for  case  4  at  the  two  different  current  densities  of  0.2  A  cm-2  and 
1.2  A  cm-2.  The  interfacial  droplet  coverage  at  the  cathode  GDL/GC,  Sint  is  assumed  to 
be  15%  for  the  lower  current  density  case  (0.2  A  cm-2)  and  7.5%  for  the  higher  current 
density  case  (1.2  A  cm-2).  The  simulation  conditions  were  at  80  °C  and  1.5  atm  for  both 
anode  and  cathode  GC  inlet  pressures. 


cathode  GDL  despite  the  significant  difference  in  operating  current 
density  and  thus  water  generation  rate  among  these  cases.  The 
numerical  trend  is  almost  identical  to  the  trend  observed  in  the  NR 
data  measured  by  Manahana  et  al.  [1  ]  at  0.2  A  cm-2  and  1.2  A  cm-2 
(see  Fig.  la). 

Our  comparisons  clearly  indicate  that  an  overall  characteristic  of 
liquid  water  transport  through  a  GDL  under  inhomogeneous 
wetting  and  porous  conditions  can  be  macroscopically  approxi¬ 
mated  by  the  proposed  two-phase  model  with  the  consideration  of 
contact  angle  dependence,  i.e.  the  second  term  on  the  right-hand 
side  of  Eq.  (18).  This  work  successfully  enhances  the  capability  of 
a  computational  two-phase  fuel  cell  model  for  a  realistic  prediction 
of  the  liquid  water  profile  and  flooding  behavior  in  a  GDL. 

4.  Conclusion  and  future  plans 

The  through-plane  water  distributions  in  GDLs  measured  by 
Manahan  et  al.  [1]  and  Turhan  et  al.  [2]  exhibit  liquid  saturation 
peaks  near  the  centers  of  the  GDLs.  In  this  study,  we  demonstrated 
that  the  experimental  trend  observed  in  [1,2]  was  successfully 
captured  by  our  macroscopic  two-phase  model  in  which  we 
included  the  proposed  term,  i.e.  the  second  term  on  the  right-hand 
side  of  Eq.  (18)  describing  two-phase  transport  through  an  inho¬ 
mogeneous  GDL  with  spatial  variation  of  wettability.  Our  model 
simulation  results  indicate  that  the  centrally  located  liquid  satura¬ 
tion  peaks  in  the  GDLs  can  be  attributed  to  relatively  fewer 
PTFE-coated  pores  in  the  inner  GDL  region  than  in  the  outer  GDL 
region  due  to  incomplete  PTFE  treatment.  Since  the  two-phase 
transport  characteristic  has  a  significant  effect  on  GDL  liquid  satu¬ 
ration  profiles  and  flooding  behavior,  efforts  are  currently  underway 
to  extend  the  present  1  -D  fuel  cell  model  to  a  three-dimensional  fuel 
cell  model  in  order  to  more  comprehensively  study  the  phenomena. 
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Nomenclature 

a  water  activity  or  effective  catalyst  area  per  unit  of  total 

volume  (m2  m-3) 

Dk  mass  diffusivity  of  species  1<  (m2  s-1) 

EW  equivalent  weight  of  dry  membrane  (kg  mol-1) 

F  Faraday  constant  (96,487  C  mol-1) 

I  current  density  (A  cm-2) 

/  diffusive  mass  flux  of  i  phase  (kg  m-2  s-1). 

J  Leverett  function 

I<  hydraulic  permeability  (m2) 

m  mass  fraction 

M  molecular  weight  (kg  mol-1) 

rid  electro-osmotic  drag  coefficient 

P  pressure  (Pa) 

Pc  capillary  pressure  (Pa) 

RH  inlet  relative  humidification 

Ru  universal  gas  coefficient  (8.314  J  mol-1  K-1) 

Sim  immobile  liquid  saturation 

Sint  interfacial  liquid  coverage 

sr  reduced  liquid  saturation 

T  temperature  (I<) 

u  fluid  velocity  and  superficial  velocity  in  a  porous  medium 

(ms-1) 
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V  volume  (m3) 

Greek  symbols: 

y  advection  correction  factor 

5l  thickness  of  component  i 

e  volume  fraction  of  gaseous  phase  in  porous  region 

6  contact  angle  (°) 

A  membrane  water  content  (molH2o  m°lso3 ) 

Xa  relative  mobility  of  phase  a 

H  viscosity  (kg  m-1  s-1) 

v  kinematic  viscosity  (m2  s-2) 

p  density  (kg  rrr3) 

pmem  dry  membrane  density  (kg  m~3) 

a  surface  tension  (N  rrr1 )  or  electronic  conductivity  (S  m-1 ) 

Superscripts: 

eff  effective  value  in  porous  region 

g  gas 

1  liquid 

mem  membrane 

Subscripts: 

GC  gas  channel 

GDL  gas  diffusion  layer 

H2  hydrogen 

mem  membrane 

N2  nitrogen 

02  oxygen 

pore  pore 

w  water 
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